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The giant arc statistics in the three year WMAP 
cosmological model 
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ABSTRACT 

We use high-resolution iV-body simulations to investigate the optical depth of giant 
arcs with length-to-width ratio larger than 7.5 and 10 in the 'standard' ACDM model 
with (j § = 0.9 and fl m ,o = 0.3 and a model based on three-year Wilkinson Microwave 
Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the 
lensing probability in the three-year WMAP model (with erg = 0.74 and f2 m .o = 0.238) 
decreases by a factor of ~ 6 compared with that in the 'standard' ACDM model. The 
effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue 
that baryons will only increase the the lensing cross-section by a moderate factor, 
~ 2. We conclude that the low central value of a% and £l nii o preferred by the WMAP 
three-year data may be too low to be compatible with observations if conventional 
assumptions of the background source population are correct. 
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1 INTRODUCTION 

In the hierarchical scenario of structure formation, struc- 
tures in the universe are assumed to have grown from tiny 
quantum fluctuations generated during an inflation period 
through gravitational instability. Inflationary models predict 
that the initial fluctuations are Gaussian and have a roughly 
scale-free power-spectrum: P(k) oc k n with n ~ 1. The am- 
plitude of the power spectrum, for which a reliable prediction 
from inflationary theory is still lacking, has to be determined 
from observations. This amplitude is usually represented by 
a quantity, erg, which is the RMS of the linear density per- 
turbations within a spherical window of radius 8/i -1 Mpc, 
where h the Hubble constant in units of 100 kms -1 Mpc - . 

The linear power spectrum can be probed with a variety 
of observations, ranging from abundances of clusters, weak 
gravitational lensing, large-scale structure in galaxy distri- 
bution, galaxy motions, Lyman-a forests, and most impor- 
tantly, the temperature fluctuations in the cosmic microwave 
background (see Spergel et al. 2006, and references therein). 
Recent results from the three-year data of the Wilkinson 
Microwave Anisotropy Probe (WMAP) alone prefers a erg 
value of 0.74±o.o5- 

Although the errorbar is still quite large, 
the median value is lower than the value, as = 0.9, adopted 
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in the 'standard' ACDM model. The lower value of as is 
favored by a number of other observations related to the 
clustering of galaxies (e.g. Jing, Mo & Borner 1998; Yang et 
al. 2004, 2005; van den Bosch, Mo & Yang 2003; Tinker et 
al. 2006; Sanchez et al. 2006), but weak lensi ng results (e.g. 
iBacon et aT]l2004l IVan Waerbeke et al-lfcooj ) and data on 
Lyman alpha forests (Viel, Haehnelt & Lewis 2006; Seljak, 
Slosar & McDonald 2006) seem to prefer a higher value. 

As the largest bound structures in the universe, clus- 
ters of galaxies provide a sensitive probe for erg, because 
their abundance in the universe depends sensitively on the 
amplitude of the density perturbations (and also the mat- 
ter content in the universe, f2 m ,o)- The abundance of clus- 
ters can be determined from X-ray fe.g.lRosati et alJl2002l: 
[ Popesso et all 120041: iReiprich fc Bohringerl 120021: iReipricbl 
l2006h and optical surveys (e.g. iGladders et al]|2006T) . How- 
ever, in order to use the observed abundance to constrain 
model of structure formation, one has to know the masses 
of these objects accurately. This is in general difficult to 
do. For example, X-ray studies usually assume hydrostatic 
equilibrium in order to derive the cluster mass, which 
may be invalid for many (such as merging) clusters (e.g. 
Gao & White 2006). Another probe of the abundance of 
clusters is giant arcs that are produced when background 
galaxies are tangentially stretched by foreground clus- 
ters (e.g. iLuppino et al] 119991 IZaritskv fc Gonzalez! 120031 : 
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Table 1. Cosmological simulation parameters. The columns are 
cosmology, box size, number and mass of dark matter particles, 
and the softening length. 



model 


box size 


A D M 




softening 




(ft-^pc) 






(h- 1 kpc) 


ACDMO 


300 


512 3 


1.67 x 10 10 


30 


WMAP3 


300 


512 3 


1.32 x 10 10 


30 



iGladders etafl l2003t ISand et alJ l2005h . A number of re- 
cent investigations made comparisons between the observed 
abundance of giant arcs with that exp ected in the 'stan- 
dard' ACDM model with a 8 = 0.9 foalal et al.| 12004 " 



iLi et al.l200llHennawi et al.l2005l:lHoresh et alJ2005ft . find- 
ing agreement between model and observation iDalal et alJ 
2004; Wambs ganss et al J2004) . In light of the change in cos- 
mological parameters preferred by the WMAP three-year 
data, it becomes important to re-evaluate how the predicted 
lensing probability changes. 

In this paper, we study the number of giant arcs ex- 
pected in the new cosmological model, using high resolu- 
tion A^-body simulations. We compare the results with those 
obtained from similar simulations of the standard ACDM 
model, and with current observations. The plan of the pa- 
per is as follows. In §2, we describe the simulations we use 
and the analysis to predict giant arcs in the Af-body simu- 
lations. Our main results are presented in §3, and we finish 
with a discussion in 54. 



2 COSMOLOGICAL MODELS, NUMERICAL 
SIMULATIONS, AND LENSING METHOD 

In this paper, we use two sets of simulations, one for the 
'standard' ACDM model with <rg = 0.9, and the other for the 
cosmological model given by the recent three-year WMAP 
data with erg = 0.74. For brevity, these two models will be 
referred to as ACDMO, and WMAP3, respectively, and the 
corresponding cosmological parameters are: 

(i) ACDMO: O m ,o = 0.3,fi A ,o = 0.7, h = 0.7, a 8 = 
0.9, n = 1; 

(ii) WMAP3: O m , = 0.238, ft A ,o = 0.762, h = 0.73, a s = 
0.74,n = 0.95. 

The assumed initial transfer function in each model was gen- 
erated with CMBFAST (Seljak & Zaldarriaga 1996). Notice 
that both erg and fi m ,o differ in these two models. 

We use a vectorized-parallel P 3 M code (Jing 
fc Suto 2002) and a PM-TREE code - GADGET2 
JSpringel. Yoshida. Whitell200 ll ISpringellliooil) to simulate 
the structure formation in these models; the details of 
the simulations are given in Table 1. Both are A'-body 
simulations which evolve Adm = 512 3 dark matter particles 
in a large cubic box with sidelength equal to 300/i _1 Mpc. 
The ACDMO simulation was performed usin g the P 3 M code 
of Jing & Suto (2002) and has been used in lLi et alJ (I2005T) 
to study the properties of giant arcs. We refer the readers to 
that paper for the detail. The WMAP3 simulation is a new 
simulation carried out with GADGET2. With their large 
volumes, these two simulations sample the corresponding 
cluster mass functions reasonably well, and they will be 



used to compare the optical depths of giant arcs in the 
ACDMO and WMAP3 models. 

For each simulation, we identify virialised dark mat- 
ter halos using the friends-of-friends method with a linking 
length equal to 0.2 times the mean particle separation. The 
halo mass M is defined as the virial mass within the virial 
radius according to the spherical collapse model (Kitayama 
& Suto 1996; Bryan & Norman 1998; Jing & Suto 2002). 
Giant arcs are produced mainly by clusters of galaxies, and 
so we focus on massive haloes with M > 10 14 h~ 1 MQ. 

For a given cluster, we calcula te the smoothe d surface 
density maps using the method of ILi et ail i2006t) . Specif- 
ically, for any line of sight, we obtain the surface density 
on a 1024 x 1024 grid covering a square of (comoving) 
sidelength of 4/i -1 Mpc centered on each cluster. The pro- 
jection depth is chosen to be 4/i _1 Mpc. Particles outside 
this cube and large-scale structures do not contribute sig- 
nificantly to the len sing cross-section (e.g. ILi et al.ll2005l : 
iHennawi et all2005t) . Our projection and smoothing method 
uses a smoothed particle hydrodynamics (SPH) kernel to 
distribute the particle mass on a 3D grid and then inte- 
grat e along the lin e of sight to obtain the surface density 
fsee ILi et al.ll20o3 for detail). In this work, the number of 
neighbors used in the SPH smoothing kernel is fixed to be 
32. Once a surface density map is obtained, we compute the 
cross-sect ion of gi a nt ar c formation following the method 
given m ILi et al] J2005I) . We consider seven source red- 
shifts, z s = 0.6, 1.0, 1.5, 2.0, 3.0, 4.0 and 7.0. The background 
sources are assumed to be elliptical with random position an- 
gles and a fixed angular surface area, S S ource = 7rx0.5" 2 . The 
axis ratio is randomly drawn between 0.5 and 1. This choice 
of axis ratio and source size is in good agre ement with the 
study of high-redshift galaxies using HST bv lFerguson et al] 
i2004l) . We generate a large number of background sources 
within a rectangle box with area Sbox- This rectangle is cho- 
sen to enclose all the high-magnification regions that can po- 
tentially form giant arcs. The number of sources generated 
is given by n sourcc = 9 Sbo*/ Source- For each source, ray- 
tracing is used to find the resulting image(s). Giant arcs are 
identified as elongated images with length-to- width (L/W) 
ratio exceeding 7.5 or 10, the usual criterion used to select 
giant arcs in observations. We calculate the total cross sec- 
tions of the top 200 most massive clusters in each simulation 
output and obtain the average cross section per unit comov- 
ing volume by: 



Zl, Z K 



V 



(1) 



where ai(z\,z s ) is the average cross-section of the three pro- 
jections of the i-th cluster at redshift z\, z s is the source 
redshift, and V is the comoving volume of the simulation 
box. The optical depth can then be calculated as: 



t(z s ) = 



4ttD s 



dz a(z, z s )(l + z) 



dV p (z) 
dz 



(2) 



where D a is the angular diameter distance to the source, and 
dVp(z) is the proper volume of a spherical shell with redshift 
from z to z + dz. The integration step size is the same as 
the redshift interval of simulation output (dz « 0.1). 
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Figure 1. The left panel shows the number density of halos in 
the ACDMO (thick) and WMAP3 (thin) cosmological models (see 
text) at redshift 0.33 obtained using the Sheth-Tormen mass func- 
tion (2002). The thin solid histogram shows the number density 
of haloes found in a simulation with box size 300ft,- 1 Mpc in the 
WMAP3 model at z = 0.33. The thick solid histogram shows the 
number density of halos in the ACDMO simulation at z = 0.2896 
(the corresponding theoretical prediction is shown as the thick 
dashed line). The right panel shows how the number density of 
haloes with mass M > 10 14 /i _1 Mq evolves as a function of red- 
shift in the two cosmologies. 



3 RESULTS 

The left panel of Fig.0sb.ows the mass functions in the two 
cosmological models at redshift 0.3, the optimal lensing red- 
shift for a source at redshift 1 (see the right panel of Fig.|2J. 
The mass functions are calculated following Sheth & Tor- 
men (2002), which is based on the Press-Schechter (1974) 
formalism and the ellipsoidal collapse model (Sheth, Mo & 
Tormen 2001). The mass functions (shown as histograms) in 
simulations clearly match the predictions well. We can see 
that the number density of haloes in the WMAP3 model is 
always smaller than that in the ACDMO model on cluster 
mass scale. For M ~ 1O 14 /i _1 M0, the abundance of haloes 
is lower by a factor of 2.6 in the WMAP3 model compared 
with that in the ACDMO model; for M ~ 10 15 /i _1 M©, the 
reduction factor is ~ 8.3. The lower abundance of clusters in 
the WMAP3 model implies that ideally a larger simulation 
box is required to sample a similar number of massive clus- 
ters as in the ACDMO model, a point we return to briefly in 
the discussion. 

The right panel in Fig.Qshows how the predicted num- 
ber density of haloes above 10 14 /i _1 Mq evolves as a function 
of redshift in these two cosmologies. The ratio of the number 
density of massive haloes in the WMAP3 model and that in 
the ACDMO model as a function of redshift is well approxi- 
mated (within 10%) by 0.4-0.2 z for z < 1.5. At redshift 0.3, 
0.5 and 1, the number density of haloes above 10 14 /i _1 Mq 
is about 1/3, 1/5 and 1/10 of that in the ACDMO respec- 
tively. Notice that the cosmological distances and volumes 
also differ in the three models, but the differences (which are 
accounted for in our calculations) are small, ranging from 
11% for the angular diameter distances and 6% for the vol- 
ume out to a source redshift of 2. These differences are much 
smaller than the change in the cluster abundance. 

As a consequence of the lower cluster abundance in 
the WMAP3 model, one expects the lensing optical depth 
will also be much lower compared with that in the ACDMO 
model. This is illustrated in the left panel of Fig. [5] which 
shows the optical depths as a function of the source redshift 
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Figure 2. The left panel shows the optical depth as a function 
of source redshift in the ACDMO and WMAP3 models for giant 
arcs with L/W > 7.5 and L/W > 10. The right panel shows 
the differential probability distribution for the optical depth as 
a function of lens redshift for z s = 1,2 and 7 respectively. In 
both panels, the results for the ACDMO and WMAP3 models are 
shown as thick and thin curves respectively. 



in the ACDMO and WMAP3 models. First, we notice that 
the optical depth is a strong function of the source redshift, 
z s , particularly when z 3 < 3, a trend first emphasized by 
IWambseanss et"al] i2004l) . Second, it is clear that the opti- 
cal depth in the WMAP3 model is much lower: compared 
with that in the ACDMO model, the optical depth is re- 
duced by more than a factor of ~ 6 for L/W > 7.5 and 
10 at all source redshifts. The right panel of Fig. [5] shows 
the differential probability distribution of the optical depth 
for giant arcs with L/W > 7.5 as a function of the lens 
redshift for z s — 1, 2, and 7; the trends (not shown) are sim- 
ilar for L/W > 10. For a source at redshift 1, 2 and 7, the 
optimal lensing redshift is around 0.3, 0.5 and 0.75, respec- 
tively, in the two cosmologies considered here. Note that the 
normalised shapes are quite similar in the two cosmological 
models but the overall optical depth is much lower in the 
WMAP3 model. 



4 DISCUSSIONS 

We have compared, using high-resolution TV-body simula- 
tions, the lensing probabilities of giant arcs in two cosmolog- 
ical models: the old 'concordance' model with as = 0.9, and 
the model with parameters given by the three-year WMAP 
data. We find that the lensing probability decreases by a 
factor of ~ 6 in the 3-year WMAP model compared to that 
in the ACDMO concordance model. This decrease is largely 
a result of the much lower number density of massive haloes 
in the WMAP3 model (see Fig. 1). 

There has been a long debate whether the number 
of observed giant arcs is consistent with theoretical pre- 
dict ions. Early comparisons use simple analytical mod- 
els (Kormann ct al ]|1994JWu fc Hammerlll99llWu fc Mai 
Il996t see also lQeuri et alJl2003l) . which were later shown to 
severely under-estimate the lensing optical depth (Bartel- 
mann & Weiss 1994; Bartelmann, Steinmetz & Weiss 1995; 
Meneghetti et al. 2003a; Torri et al. 2004). Many re- 
cent studies use N-body simulations, but with only d ark 
matter (e.g. IWambseanss et alJl20o5 : iDalal et alJl2004l . see 
iPuchwein et al.l2003 for an exception) . Below we discuss the 
uncertainties in both observations and predictions, and ex- 
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amine the consistency found in several previous works in the 
ACDMO cosmology. 

4.1 Consistency in the ACDMO cosmology 

Several previous studies concluded that the observa- 
tional results o f giant arcs are compatible with theoreti- 
cal predictions fOeuri et alJl2003t IWambseanss et aLll2004i 
iDalal et alJl2004) in the ACDMO cosmology. However, there 
are uncertainties in such conclu s ion. 

The study by lOeuri et alJ J2003I) used the axis ratio 
distributions in the tri-axial model of Jing & Suto (2002), 
and assu med a central profile, p oc r~ 1 5 , that are steeper 
than the lNavarro. Frenk fc White! il997t) form. Gas cooling 
can steepen the halo density profile, but it also makes the 
central mass distribution rounder. In this case, the shape 
distributions obtained in dark matter-only simulations are 
not suitable for clusters where gas cooling steepens the den- 
sity profile dramatically. Since a shallower profile or a more 
spherical shape will reduce the giant arc cros s-sections, it 
is unc lear whether or not the conclusion of lOeuri et alJ 

(2003) still holds if con sistent models for densi ty profile 
and halo shape are used. Wambsaanss ct al. (2004) assumes 
that the L/W ratio can be approximated by the magni- 
fication. This assumption is valid for isothermal spheres, 
but not for real clusters which appear t o be well fit by the 

[ Navarro. Frenk fc White! dl997l) profile (Ivan de Marel et alJ 
20001: IComerford et aljl2006l: IVoiet fc Fabianl2006l) . This as- 
sumption over-estimates the optic al depth by a factor of few 
JPalal et alJ 12004 iLi et aljl2005l) . Furthermore, they also 
assumed a slightly higher normalization, ag( 0.95), which 
also in creases the optical depth. The study bv lDalal et alJ 

(2004) assumes that arcs can be approximated as rectangles, 
which may overestimate the total cross-section by a factor 
of two compared to the perhaps more realistic assumption 
of elliptical arcs. In summary, even in the ACDMO cosmol- 
ogy, the consistency between observations and predictions 
requires somewhat optimistic assumptions. The WMAP3 
model, which predicts an optical depth by a factor of ~ 6 
smaller, makes it even harder to explain the observed giant 
arcs. 



4.2 Effects of baryons 

In the real universe, baryons account for roughly 20% of 
the total mass, which can cool (form stars) and sink toward 
the centres of clusters. The radiative cooling likely has two 
effects: it will increase the concentration of baryons at the 
center of clusters, and at the same time, make the cluster s 
more spherical ('e.g. lDubinskil9'94llKazantzidis et all2004Tl . 
The former increases while the latter decreases the lensing 
cross-sections, and so the overall influence depends on which 
effect dominates. 

Puchwein et al. (2005) studied the lensing cross-sections 
of clusters with M > 10 15 Mq with different prescriptions 
of numerical viscousities. They concluded that the higher 
concentrations due to baryons dominate and the baryons in- 
crease the lensing cross-sections by a factor of ~ 2. However, 
they cautioned that their clusters may suffer from overcool- 
ing, as the stellar density in the core is larger than observed, 
and so the effect of baryons may be over-estimated. Notice 



also that their simulations did not include feedbacks from ac- 
tive galactic nuclei, which may further decrease the cooling 
of baryons at cluster centres and the corresponding lensing 
cross-sections. As the star format ion treatment in hydrody - 
namical simulations is uncertain. iMeneehetti et al] (|2003b) 
adopted an alternate empirical approach. They studied the 
effect of baryons by including a central cD galaxy in clusters. 
They found that the increase in the lensing cross-section due 
to cD galaxies is again quite moderate, by a factor of < 2. 
Such an enhancement cannot compensate the reduction in 
the optical depth due to the lower abundance of clusters in 
the WMAP3 model compared with the ACDMO model. 

4.3 Uncertainties in the cluster and background 
source populations 

The cluster abundance (but not internal structures, partic- 
ularly when star formation is included) can now be reliably 
predicted from both numerical simulations and the extended 
Press-Schechter analytic formalism. In the ACDMO cosmol- 
ogy, various investigations, using similar assumptions of the 
background source population and definition of the length- 
to- width ratio, predict lensing optical dep ths that agree with 
one another w ithi n a factor of 1.5 (e.g. IDalal et al] 12003 : 
Li et al. 2005; sec ILi et aljl2005l for a detailed discussion). 
We caution, however, that even our 300ft" 1 Mpc box size 
still appears to be somewhat too small to sample the tail of 
the cluster mass function well, as is shown by the large noise 
in the right panel of Fig. 2. We have 158 clusters with M > 
10 14 /i _1 Me (7 clusters with M > 3.8 x lO^/i^M©). At red- 
shift z = 0.3 for z s — 3, 50% of the giant arc cross-section 
is contributed by clusters with M > 3.8 x 1O 14 /i _1 M . The 
larger noise may also be partly due to our inadequate sam- 
pling of merger events as we only dump the simulation data 
with a redshift interval of dz ~ 0.1. As shown by a num- 
ber of previous works (Torri et al 2004; Hennawi et al. 2005; 
Fedeli et al. 2006), merger events can boost the lensing cross- 
sections significantly. To examine the importance of this ef- 
fect, we ran a lower resolution simulation in the WMAP3 
model which evolves A?dm = 512^ dark matter particles in 
a box with sidelength of 600/i~ Mpc. Due to the larger vol- 
ume, the new simulation better samples the high mass tail 
of the cluster mass function and merger events. The optical 
depth decreases by a factor of ~ 2. The larger optical depth 
in our first WMAP3 simulation is partly because it has rel- 
atively more massive clusters than the larger simulation due 
to cosmic variance. Furthermore, the poorer resolution of the 
larger simulation may have also reduced the optical depth, 
and so the real change should be < 2. 

Substructures along the line of sigh t may be impo rtant 
for the anomalous flux ratio problem dMetcalJ 12005). but 
are unlikely to be important for the lensing cross-sections 
as they are expected to contribute only a few per cent of 
the surface density. Notice that we use a cube of 4ft Mpc 
to evaluate the lensing cross-sections. In lLi et ail (12005!) the 
sidelength is chosen to be 2r v i r . These two choices give al- 
most identica l optical depths (compare Fig. 2 with Fig. 7 in 
ILi et al.l2005T) . This demonstrates that the influence of mat- 
ter (including substructures) in the outer skirts of clusters 
is not important. 

A much larger uncertainty concerns the background 
source population, including their ellipticity, size and red- 
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shift distributions. While the former two app ear to have 
modest effect s on the lensing cross-sections fe.g. lOguril2002l : 
iLi et, al.l2005F ) , the source redshift distribution has quite dra- 
matic effects on the lensing cross-section (Wambsganss et al. 
2004), as can be seen in Fig. [5] This is the biggest uncer- 
tainty in the predictions of giant arcs. More redshift mea- 
surements of gian t arcs will be partic ularly useful to clarify 
the situation fsee lCovone et al . 2006 for a recent effort). 

Currently the giant arc samples are still small. The 
largest dedicated search for giant arcs in X-ray clusters was 
performed by Luppino et al. (1999) who found strong lensing 
in 8 out of 38 clusters. In the optical, Zaritsky & Gonzalez 
(2003) found two giant arcs using the Las Campanas Distant 
Cluster Survey between redshift 0.5 and 0.7 in an effective 
search area of 69 square degrees while Gladders et al. (2003) 
found eight clusters with giant arcs using the Red-Sequence 
Cluster Survey over ~ 90 square degrees. The largest arc 
sample found from HST archives was presented by Sand et 
al. (2005) with 116 arcs from 128 clusters, although its selec- 
tion function is likely heterogeneous. Clearly the current gi- 
ant arc samples are somewhat limited, alt hough future weak 
lensing surveys (e.g. IWittman et al JI2006I) may yield a large 
number of giant arcs as a by-product. A combined analy- 
sis of strong and weak lensing will be particularly useful for 
constraining as and the cluster inner mass profiles. 

Until larger arc samples become available and our un- 
derstanding of the background source population is im- 
proved, it is difficult to reach firm conclusions concerning 
the consistency of the WMAP three-year model with ob- 
servations. Nevertheless, it appears difficult to reconcile the 
giant arc statistics with the low central <T8 value (0.74) pre- 
ferred by the WMAP three-year data. 
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